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ABSTRACT 

We use a hybrid, multiannulus, n-body-coagulation code to investigate the 
growth of km-sized planetesimals at 0.4-2 AU around a solar-type star. After 
a short runaway growth phase, protoplanets with masses of ~ 10^® g and larger 
form throughout the grid. When (i) the mass in these ‘oligarchs’ is roughly 
comparable to the mass in planetesimals and (ii) the surface density in oligarchs 
exceeds 2-3 g cm“^ at 1 AU, strong dynamical interactions among oligarchs 
produce a high merger rate which leads to the formation of several terrestrial 
planets. In disks with lower surface density, milder interactions produce several 
lower mass planets. In all disks, the planet formation timescale is ~ 10-100 Myr, 
similar to estimates derived from the cratering record and radiometric data. 

Subject headings: planetary systems - solar system: formation - stars: formation 
- circumstellar matter 


1. INTRODUCTION 

In the protosolar nebula, the growth of terrestrial planets begins with collisions and 
mergers of planetesimals, solid objects with radii of ~ 1 km (e.g. Safronov 1969; Lissauer 


^Calculations reported here used the ‘hydra’ cluster run by the Computation Facility of the Harvard- 
Smithsonian Center for Astrophysics 
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1987; Wetherill & Stewart 1993; Weidenschilling et al. 1997). Initially, planetesimals 
grow slowly. As they grow, dynamical friction circularizes the orbits of the largest objects; 
viscous stirring excites the orbits of the smallest objects. Slow, orderly growth ends when 
the gravitational cross-sections of the largest objects exceed their geometric cross-sections. 
Because dynamical friction is faster than accretion, the largest objects stay on circular orbits 
and grow faster and faster relative to the smallest planetesimals. This runaway growth 
rapidly concentrates solid material into a few protoplanets (e.g., Greenberg et al 1978; 
Greenberg et al. 1984; Wetherill & Stewart 1989, 1993; Weidenschilling et al. 1997; 
Kenyon & Luu 1998; Kokubo & Ida 2000). 

During runaway growth, protoplanets stir the orbits of the leftover planetesimals. Stir¬ 
ring reduces gravitational focusing factors and slows the growth of the protoplanets (Wetherill 
& Stewart 1993; Ida & Makino 1993; Kenyon & Bromley 2002; Rahkov 2003c). Although 
protoplanets grow slowly, they grow faster than the leftover planetesimals and intermedi¬ 
ate mass objects do. These large ‘oligarchs’ slowly clear their orbits of smaller objects and 
reach a maximum ‘isolation mass’ that depends on the initial surface density of solid ma¬ 
terial (Lissauer 1987; Lissauer & Stewart 1993; Kokubo & Ida 1998; Rahkov 2003c; 
Goldreich, Lithwick & Sari 2004b). 

Throughout oligarchic growth, protoplanets repel other oligarchs in the disk (Kokubo 
& Ida 1996, 1998, 2000; Kominami & Ida 2002; Thommes, Duncan, & Levison 2003). 
Eventually, however, dynamical interactions produce collisions and gravitational scatter¬ 
ing among the oligarchs (Ghambers 2001; Kokubo & Ida 2002; Kominami & Ida 2002; 
Thommes, Duncan, & Levison 2002; Rahkov 2003c; Goldreich, Lithwick & Sari 2004b). 
During this phase of ‘chaotic growth’, the largest oligarchs merge with and clear away other 
oligarchs and smaller objects to become full-hedged planets (e.g., Goldreich, Lithwick & Sari 
2004a). 

Within this framework, the transition from oligarchic growth to chaotic growth is poorly 
understood. Although analytic estimates provide a guide to the late evolutionary stages (e.g., 
Rahkov 2003c; Goldreich, Lithwick & Sari 2004a), numerical calculations are necessary to 
test the basic theory and to derive the end states and timescales as a function of initial 
conditions (e.g. Kominami & Ida 2002). Published calculations do not test this regime of 
the theory in much detail. Pure coagulation and simple hybrid calculations cannot follow 
the transition accurately (Wetherill & Stewart 1993; Weidenschilling et al. 1997). Direct 
n-body calculations with planetesimals are computationally expensive and tend to focus 
on evolution during oligarchic growth, when orbits of individual objects are easier to track 
(Alexander & Agnor 1998; Kokubo & Ida 2002, and references therein). Direct n-body 
calculations without planetesimals follow dynamical interactions after the transition and 
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concentrate on the clearing phase (Chambers 2001). Thus, understanding this evolutionary 
phase requires new calculations (see the discussion in Kokubo & Ida 2002). 

Here, we use numerical calculations with a new hybrid n-body-coagulation code to 
investigate the transition from oligarchic growth to chaotic growth. Our approach allows 
us to combine statistical algorithms for the planetesimals with direct n-body integrations of 
the oligarchs. From several simple test cases and complete planet formation calculations, we 
show that oligarchic growth becomes chaotic when the orbits of oligarchs begin to overlap. 
If the surface density in oligarchs exceeds a critical value, this transition occurs when the 
oligarchs contain roughly half of the mass in solids. At 1 AU, the critical initial surface 
density is T^c ~ 2-3 g cm“^. Thus, oligarchs can make the transition from oligarchic to 
chaotic growth in disks with masses comparable to the minimum mass solar nebula, where 
S ^ 8-10 g cm“^ at 1 AU. 

In disks where the surface density of solids is below the limit for chaotic growth, oligarchs 
slowly merge to form larger objects. Pairwise interactions, instead of large-scale chaos, drive 
the dynamics of these systems. Milder, slower interactions between oligarchs then produce 
less massive planets. 

We develop the background for our calculations in §2, describe a suite of calculations in 
§3, and conclude with a brief summary and conclusions in §4. 


2. BACKGROUND 

The evolution from planetesimals to planets is marked by several phases - orderly 
growth, runaway growth, oligarchic growth, and chaotic growth - with clear transitions 
in the dynamics and mutual interactions of the ensemble of solid objects. Analytic deriva¬ 
tions and sophisticated coagulation and n-body calculations identify the physics of these 
transitions. Here, we summarize some basic results to provide the context for our numerical 
simulations. 

Most considerations of planet formation begin with small objects, Tj < 1 km, that 
contain all of the solid material. For these sizes, collisional damping and viscous stirring 
roughly balance for orbital eccentricity e ~ 10“^. The gravitational binding energy. Eg ~ 10^ 
erg g“^, is then comparable to the typical collision energy at 1 AU, Ec ~ 10^-10^ erg g“^. 
Both energies are smaller than the disruption energy ~ the collision energy needed to remove 
half of the mass from the colliding pair of objects - which is Qd ^ 10^ erg g“^ for rocky 
material (Davis et ah 1985; Benz & Asphaug 1999). Thus, collisions produce mergers 
instead of debris. 
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Initially, growth from mergers is slow. The collision cross-section is the product of the 
geometric cross-section and the gravitational focusing factor, 

dc ~ Tir'^fg ~ 7rr^(l -F /3{evK/vescf) , (1) 

where r is the particle radius, vk is the orbital velocity, Vesc is the escape velocity, and 
/5 is a coefficient that accounts for three-dimensional orbits in a rotating disk (Greenzweig 
& Lissauer 1990; Spaute et ah 1991; Wetherill & Stewart 1993). Because evx ~ Vesc, 
gravitational focusing factors are small. Thus, growth is slow and orderly (Safronov 1969). 

As larger objects form, the smaller objects effectively damp the orbital eccentricity of 
larger particles through dynamical friction (e.g. Wetherill & Stewart 1993; Kokubo & Ida 
1995; Kenyon & Luu 1998). Viscous stirring by the large objects heats up the orbits of the 
small objects. In the case where gas drag is negligible, Goldreich, Lithwick & Sari (2004b) 
derive a simple relation for the ratio of the eccentricities of the large (‘1’) and the small (‘s’) 
objects in terms of their surface densities E; ^ (see also Kokubo & Ida 2002; Rahkov 2003a, 
2003b, 2003c), 



with n = 1/4 to 1/2. For S;/Ss ~ 10 ^ to 10 eijcs ~ 0.1-0.25. Because CsVk ^ 
gravitational focusing factors are large. Runaway growth begins. 

Runaway growth concentrates more and more mass in the largest objects. Dynamical 
friction produces the largest gravitational focusing factors among the largest objects. These 
protoplanets run away from the smaller objects and contain an ever growing fraction of the 
total mass. At the same time, the large objects continue to stir the leftover planetesimals. 
The leftovers have orbital velocity dispersions comparable to the escape velocities of the oli¬ 
garchs. With CsVk ~ Vesci equations (1) and (2) show that collision rates decline as runaway 
growth continues. The ensemble of protoplanets and leftover planetesimals then enters the 
oligarchic phase, where the largest objects grow faster than the leftover planetesimals. 

Among the oligarchs, the smaller oligarchs grow the fastest. Each oligarch tries to 
accrete all of the material in an annular ‘feeding zone’ set by balancing the gravity of neigh¬ 
boring oligarchs. Within each annulus, each oligarch stirs up the remaining planetesimals 
along its orbit. Because smaller oligarchs orbit in regions with smaller E^/S^, equations (1) 
and (2) show that smaller oligarchs have larger gravitational focusing factors. Thus smaller 
oligarchs grow faster (e.g., Kokubo & Ida 1998; Goldreich, Lithwick & Sari 2004b). 

During oligarchic growth, protoplanets become isolated from their surroundings (Lis¬ 
sauer & Stewart 1993; Kokubo & Ida 1998, 2000; Thommes, Duncan, & Levison 2003; 
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Goldreich, Lithwick & Sari 2004a,b). As oligarchs accrete smaller planetesimals, dynami¬ 
cal interactions push them apart to maintain a typical separation of 5-10 Hill radii, Rh = 
where m is the mass of an oligarch and ttiq is the mass of the Sun. The sep¬ 
aration depends weakly on the mass and semimajor axis of the protoplanet and the local 
surface density (Kokubo & Ida 1998, 2000). 

Oligarchic growth limits the mass of the largest protoplanets (Lissauer 1987; Lissauer 
& Stewart 1993; Kokubo & Ida 1998). If an oligarch accretes all planetesimals within a 
torus of width aacc = nRn, it has a mass m ~ 0.005-0.01 where m 0 is the mass 

of the Earth. For n ~ 5, the so-called isolation mass is niiso ~ 0.1 m 0 . Unless oligarchs 
move out of their feeding zones or additional material is brought into the feeding zone, they 
cannot grow beyond the isolation mass. 

A transition from oligarchic growth to chaotic growth appears to provide the best way 
for oligarchs to evolve into planets (e.g., Kominami & Ida 2002; Goldreich, Lithwick & 
Sari 2004b). Small oligarchs on roughly circular orbits prevent radial drift by repelling 
planetesimals and opening up a gap in the feeding zone (Rahkov 2001, 2003a,b,c,d). Larger 
oligarchs that stir up planetesimals outside the gap can grow up to the isolation mass, 
but no further (Rahkov 2003c; Goldreich, Lithwick & Sari 2004a,b). Once the mass in 
the oligarchs exceeds the mass in leftover planetesimals, dynamical interactions between 
the oligarchs move them outside their feeding zones. The oligarchs then compete for the 
remaining planetesimals (Goldreich, Lithwick & Sari 2004a,b). Goldreich, Lithwick & Sari 
(2004a) propose that large-scale dynamical interactions begin when S; > S*. However, this 
regime has not been investigated in detail by numerical simulations. 

To evaluate the conditions necessary for chaotic orbits to produce terrestrial planets, 
we calculate the formation of planets from 1-10 km planetesimals. We start with several 
test calculations to demonstrate and to examine the physical processes involved in oligarchic 
growth and the transition from oligarchy to chaos. We then consider the growth of planets 
in a small torus at 0.84-1.16 AU and a large torus at 0.4-2 AU. 

To provide a measure of the orbital interactions in the calculations, we dehne two 
parameters for the ensemble of oligarchs. From published n-body calculations, we expect 
signihcant orbital interactions between oligarchs when their typical separations are less than 
5-10 Rh (Kokubo & Ida 1995, 1998; Thommes, Duncan, & Levison 2003; Bromley & 
Kenyon 2006). With N equal mass oligarchs spaced evenly in radial distance within a torus 
of width Aa, oligarchs should begin to interact when 5 — IQNRh ~ Aa. Generalizing this 
idea to N oligarchs of any mass, we expect signihcant orbital interactions among oligarchs 
when the sum of their Hill radii is ~ 0.1-0.2 Aa. To provide a measure of the onset of 
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interactions, we define a Hill parameter as the normalized sum of Hill radii for all oligarchs: 


Ph 



^out (^in 


( 3 ) 


where RH,i is the Hill radius for each of N oligarchs and {aout) is the inner (outer) radius 
of the coagulation grid. We expect signihcant dynamical interactions due to orbit overlap 
when Ph ^ 0.1. 

The Hill parameter is useful, because we can relate pn to the surface density required 
for overlapping orbits. If rrit is the typical mass of an oligarch in equation (3), the surface 
density in oligarchs is roughly 


S; « 2.5 



m 

1026 g 


2/3 

g cm 


( 4 ) 


If oligarchic growth ends when pu ~ 0.1, equation (3) yields E; ~ 2-3 g cm which is 25% 
to 33% of the mass in a minimum mass solar nebula (Weidenschilling 1977; Hayashi 1981). 


To provide a second measure of the transition from oligarchic to chaotic growth, we 
derive an orbit crossing parameter. We dehne the absolute value of the difference in the 
semimajor axis of two oligarchs, asep,ij = l^i — For each oligarch with j % i, we evaluate 
Xij = {dsep — cijCj)/Rn^ij, where Cj is the orbital eccentricity and 


f mi + rrij \ f ai + aj \ 

(—j 


( 5 ) 


is the mutual Hill radius. We hnd the minimum value of Xij, Xmin,i = rnin{xij). The orbit 
crossing parameter is then 

/ V?'—1 


Po = 


mtot 


( 6 ) 


where mtot is the total mass in oligarchs. When po > 5, the orbits of oligarchs do not cross. 


When Po ~ 0, most orbits cross. We expect signihcant orbital interactions, mergers, and 
chaotic growth when Po ^ 0. 


3. NUMERICAL CALCULATIONS 
3.1. Methods 

To explore the end of oligarchic growth, we consider numerical calculations of planet 
formation. Because the initial number of planetesimals is large, ~ 10® to 10®, a statistical 
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approach is the only plausible method to derive the collisional growth of planetesimals (e.g. 
Safronov 1969; Wetherill & Stewart 1993). When a few objects contain most of the mass, 
the statistical approach fails. N-body methods can then treat the behavior of the largest 
objects but cannot follow the evolution of the leftover planetesimals. Here, we use a hybrid 
n-body-coagulation code, which combines the statistical and n-body approaches to follow 
the growth of 1~10 km planetesimals into planets. 

Kenyon & Bromley (2004a, and references therein) summarize our multi-annulus co¬ 
agulation code. The model grid contains N concentric annuli with widths Sai centered at 
heliocentric distances a^. Each annulus contains n{mik,t) objects of mass rriik with orbital 
eccentricity eik{t) and inclination iik{t) in M mass batches. Our solution of the coagulation 
and Fokker-Planck equations treats elastic and inelastic collisions between all particles in all 
annuli. We adopt collision rates from kinetic theory - the particle-in-a-box method - and 
use an energy-scaling algorithm to assign collision outcomes. We derive changes in orbital 
parameters from gas drag, dynamical friction, and viscous stirring (Adachi et al. 1976; 
Ohtsuki, Stewart, & Ida 2002). 

Bromley & Kenyon (2006) describe the n-body code and our methods for combining 
n-body calculations with the coagulation code. When an object in the coagulation code 
reaches a preset mass, it is ‘promoted’ into the n-body code. To set the initial orbit for this 
object, we use the three coagulation coordinates, a, e, and f, and select random values for 
the longitude of periastron and the argument of perihelion. Because the annuli have hnite 
width 5aj, where i is the index, we set the semimajor axis of the promoted object to = 
Qi -|- (0.5 — x)Sai, where a; is a random number between 0 and 1. When two or more objects 
within an annulus are promoted to the n-body code during the same timestep, we restrict the 
choices of the orbital elements to minimize orbital interactions between the newly promoted 
n-bodies. 

We follow the mutual interactions of the ensemble of ‘n-bodies’, including mergers, using 
a robust set of integrators. These integrators include accretion and drag terms that link the 
evolution of the planetesimals to the evolution of the n-bodies. To compute how rapidly 
n-bodies accrete planetesimals, we adopt a particle-in-a-box formalism. Direct calculations 
for gas drag and Fokker-Planck gravitational interactions provide rates for changes in a, e, 
and i for each n-body due to interactions with planetesimals. 

Bromley & Kenyon (2006) describe tests of the hybrid code along with initial results for 
terrestrial planet formation. When the Hill radius of a promoted n-body is small compared 
to the separation of annuli in the coagulation grid, the hybrid code reproduces published 
results of previous investigations (e.g. Greenzweig & Lissauer 1990; Wetherill & Stewart 
1993; Weidenschilling et al. 1997; Duncan, Levison & Lee 1998; Chambers 2001). In the 



calculations described below, we set the promotion mass to rupro = 1-30 g, which 

yields a small Hill radius, ~ 0.0005-0.002 a, compared to the separation of adjacent annuli 
in the coagulation grid, ~ 0.01-0.04 a. 

Our numerical calculations begin with a mass distribution of planetesimals in 32-64 
concentric annuli with initial surface density S = Eq (a/1 AU)"^, with n = 1-1.5. The spacing 
between successive mass batches is 5 = mj+i/m*. We adopt 5 = 1.4-2.0 for most calculations. 
All planetesimals start with initial eccentricity Cq and inclination As planetesimals evolve, 
the average mass rriik and orbital parameters 6*^ and iik of each mass batch change. We add 
mass batches as planetesimals grow in mass and reserve 8-16 mass batches in each annulus 
for n-bodies. Throughout the calculation, gas drag and collisions transport planetesimals in 
the coagulation code from one annulus to another. In addition to these processes, mutual 
gravitational interactions can scatter n-bodies into different annuli. 

For these calculations, we assume a simple treatment of the gas. The midplane density 

is 

Pg{a,t) = po{a)e-^/^\ (7) 

where po{a) is the initial gas density and to is a constant. We adopt a gas surface density 
Eg = lOOE and set po = Sg/277, where H is the gas scale height (Kenyon & Hartmann 
1987). Consistent with observations of pre-main sequence stars, we adopt to = 1 Myr (see, 
for example, Haisch, Lada, & Lada 2001; Young et al 2004; Calvet et ah 2005; D’Alessio 
et al. 2005). 


3.2. Simple model 

To isolate the important physical parameters in the evolution of oligarchs, we first 
consider an ensemble of equal mass planetesimals in 32 concentric annuli at 0.84-1.16 AU. 
The planetesimals have mass tUs, surface density E = Es(a/1 AU)“^/^, initial eccentricity Cq 
= 10“® and inclination ig = eo/2. Within this swarm, we embed N oligarchs with mass mi 
and surface density E = E;(a/1 AU)“^/^. The oligarchs have the same initial eccentricity 
and inclination as the planetesimals. 

To evolve this ensemble in time, we calculate gravitational stirring for all interactions. 
We allow oligarchs - but not planetesimals - to collide and merge. This assumption allows us 
to focus on the dynamical evolution of the oligarchs in the absence of planetesimal accretion. 

Figure 1 shows the evolution of the eccentricity for three sets of oligarchs with E^/E^ = 
0.5 and E; = 1-4 g cm“^. The planetesimals have m* = 10^^ g; the oligarchs have mi = 10^® g, 
comparable to the isolation mass for this grid. The Hill parameter increases from pjj = 0.04 
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(lower panel) to pn = 0.07 (middle panel) to ph = 0.13 (npper panel). In each frame, the 
colored tracks indicate the eccentricities of oligarchs in the grid. Along the track of a single 
oligarch, the color changes when two oligarchs merge. Although a single color does not track 
the motion of an individual oligarch throughout the evolution, the ensemble of curves tracks 
the merger history of the hnal set of oligarchs. Tracking backwards in time along connected 
curves yields the evolution of one of the oligarchs that remains at the end of the calculation. 
The legends list S; and the initial and hnal number of oligarchs (Figure 2). 

The eccentricity evolution is sensitive to the initial mass in the swarm. For S; = 1 
g cm“^, the oligarchs have a typical separation of 15-20 Rh and do not have signihcant 
interactions. As the oligarchs stir up the planetesimals, dynamical friction maintains a 
constant ratio e^/c; ~ (m^/m^)^/^ ~ 0.1. Although the orbits of the oligarchs become more 
and more eccentric, the growth of ei is slow (Figure 1, lower panel). It takes ~ 1 Myr to 
reach ei ~ 0.01 and another 8-9 Myr to reach e; ~ 0.02. After ~ 100 Myr, when e ~ 0.04, 
the orbits begin to overlap. 

As S/ increases, orbital interactions occur on shorter timescales (Figure 1, middle and 
upper panels). For = 2 g cm“^, the initial separation of the oligarchs is ~ 10-12 Rh- 
It takes only 5 x 10"^ yr for the typical eccentricity to reach e ~ 0.01, when the minimum 
separation between oligarchs is only 8 Rh- At ~ 10® yr, 0.02 and orbits overlap. Two 
oligarchs merge at ~ 2 x 10® yr; two more merge at ~ 1 Myr. When we stop the calculation 
at 1 Myr, two oligarchs remain in eccentric orbits, e > 0.025 and are likely to merge with 
other oligarchs. 

For S; = 4 g cm“^, orbits begin to overlap in ~ 10"^ yr. After a single merger at ~ 10"^ 
yr, orbits continue to grow more eccentric. At 5 x 10"^ yr, all orbits overlap and the merger 
rate accelerates. There are two additional mergers at ~ 10® yr, another two by 3 x 10® 
yr, and three more by ~ 1 Myr. At 1 Myr, all of the remaining oligarchs have eccentric, 
overlapping orbits and many are likely to merge over the next few Myr. 

Figure 2 illustrates the chaotic behavior of the semimajor axis in these test cases. For 
S; = 1 g cm“^, the oligarchs have a constant semimajor axis for almost 100 Myr. When the 
total mass is a factor of two larger (S; = 2 g cm“^), the semimajor axes are constant for 
~ 10® yr. Once the orbits start to overlap, several oligarchs show considerable excursions 
in semimajor axis of 0.1 AU or more, ~ 30% to 40% of the grid. Two of these oligarchs 
merge with other oligarchs. For S; = 4 g cm“^, the orbits are very chaotic, with larger radial 
excursions and many mergers. 

Figures 3 and 4 show the evolution of the number of oligarchs No and the orbit crossing 
parameter po for the three cases in Figures 1-2 and a fourth case with S; = 8 g cm“^. 
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Mergers are prominent in all but the lower surface density test. In the tests where mergers 
are important, Po ^ 0 when orbits overlap and mergers begin. 

These tests conhrm our general expectations. Oligarchs stir up planetesimals along 
their orbits. Dynamical friction maintains a fixed ratio of eijcg (equation 2). Because 
the dynamical friction timescale depends on the surface density, the orbits of oligarchs in 
more massive disks overlap on shorter timescales than the orbits of oligarchs in less massive 
disks. More massive disks contain more oligarchs, leading to more chaotic orbits on shorter 
timescales. 

To provide additional tests, we consider two more sets of calculations. Models with 
hxed S; and variable S* gauge the importance of damping by leftover planetesimals. We 
hrst consider models with = 10^^ g and then examine calculations with = 10^® g. 

Figure 5 shows the evolution of semimajor axis for three calculations with = 2 g cm“^ 
(upper panel), = 4 g cm“^ (middle panel), and = 8 g cm“^ (lower panel). The tracks 
in the upper panel repeat those from the middle panel of Figure 2; chaos begins at 10^ yr 
when the orbits start to overlap. As increases, dynamical friction between the small and 
large objects is more efficient. The orbits overlap and chaos begins later when < 0.5. 

The middle set of tracks in Figure 5 demonstrates this behavior: there is no chaos until the 
system has evolved for 1 Myr. However, in models where exceeds ~ 8 g cm“^, viscous 
stirring among the planetesimals increases their velocity dispersions on shorter timescales. 
Larger viscous stirring results in larger orbital eccentricities for the oligarchs (equation 2) 
and a faster onset of chaos. For = 8 g cm“^, orbits begin to overlap at ~ 3 x 10® yr. 
These chaotic interactions grow with 

To test the importance of viscous stirring among planetesimals. Figure 6 repeats the 
calculations of Figure 5 for planetesimals with rris = 10^® g. In this test, the ratio of 
orbital eccentricity is much larger, Csjei ~ 75-80. The viscous stirring timescale for the 
planetesimals is much larger than 1 Myr. For all = 4-32 g cm“^, planetesimal stirring is 
negligible. Because dynamical friction is important, the oligarchs remain well-separated and 
never develop overlapping orbits. 

These tests demonstrate the main physical processes involved in the transition from 
oligarchy to chaos. Dynamical friction maintains a hxed ratio of eijcs (equation 2). Viscous 
stirring increases the orbital eccentricities of planetesimals until the orbits of oligarchs inter¬ 
act. When large planetesimals contain a signihcant amount of mass, they contribute to the 
stirring. Otherwise, oligarchs provide all of the stirring. Once orbits overlap, chaos ensues. 
Chaos produces mergers and large excursions in semimajor axis, which starts the process 
that clears out the disk to produce large planets. 
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Three parameters - Sz/S^, pn, and po - provide good measures of the transition from 
oligarchy to chaos. The Hill parameter measures when the oligarchs have enough mass to 
interact dynamically. The ratio isolates the time when planetesimals cannot damp the 

oligarchs and thus prevent large-scale dynamical interactions. The orbit overlap parameter 
distinguishes times when orbit overlap is important. 

To understand the transition from oligarchy to chaos in less idealized situations, we now 
consider complete planet formation simulations using the full hybrid code. The calculations 
start with 1~10 km planetesimals and allow all objects to collide, merge, and interact grav¬ 
itationally. When objects in the coagulation code reach m ~ 2 x 10^® (So/8 g cm“^) g, we 
promote them into the n-body code and follow their individual trajectories. We describe 
calculations in a small (large) torus in §3.3 (§3.4). 


3.3. Planet formation at 0.86—1.14 AU 

The calculations begin with 1-3 km planetesimals in a torus extending from 0.86 AU to 
1.14 AU. We divide this region into 32 annuli and seed each annulus with planetesimals in 
nearly circular and coplanar orbits (cq = 10“^ and io = eo/2). The planetesimals have surface 
density S = So(a/l AU)“^/^, with Sq = 1-16 g cm“^ at 1 AU. In these calculations, we do 
not consider fragmentation, which generally speeds up the growth of the largest objects at 
the expense of mass loss from disruptions and gas drag (Wetherill & Stewart 1993; Kenyon 
& Luu 1998). Weidenschilling et ah (1997) consider a similar suite of calculations. Where 
it is possible to compare, our results agree with these calculations (see also Kominami & Ida 
2002 ). 

For So = 8 g cm“^, growth at 1 AU follows a standard pattern (Wetherill & Stewart 
1993; Weidenschilling et al. 1997; Kenyon & Luu 1998). After a few thousand years, mergers 
produce a few large objects with radii of ~ 10 km. As dynamical friction circularizes the 
orbits of these objects, runaway growth begins. It takes only 10^ yr to produce several dozen 
100-300 km objects. At ~ 2 x 10^ yr, the first object is promoted into the n-body code. As 
larger objects form farther out in the disk, more promotions occur. These objects continue 
to grow rapidly until they reach ‘isolation’ masses of ~ 10^® g, when stirring begins to reduce 
gravitational focusing factors. 

The transition to oligarchic growth begins at the inner edge of the grid and rapidly 
propagates outwards. At ~ 3 x 10^ yr, the number of oligarchs with masses m > 10^® g 
peaks at No ~ 7. Soon after oligarchic growth begins at the outer edge of the grid, oligarchs 
at the inner edge of the grid begin to interact dynamically. A wave of strong dynamical 
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interactions then moves out through the grid. It takes ~ 1 Myr for the wave to move from ~ 
0.85 AU to ~ 1.15 AU. During this period, some oligarchs merge. Others migrate through 
the grid on highly eccentricity orbits. 

From ~ 1 Myr onward, mergers slowly reduce No. It takes ~ 1 Myr for the hrst 2 
mergers and another ~ 2 Myr for the second 2 mergers. After 100 Myr, only 3 oligarchs 
remain. One of these has m ~ 0.43 m®, a ~ 0.9 AU, and e ~ 0.08. The other two oligarchs 
have m ~ 0.05 uiq and e ~ 0.1 (Figure 7). Aside from the eccentricity of the more massive 
planet, the properties of these objects are reasonably close to those of the Earth and Mars. 
Fragmentation and interactions with the gas probably promote smaller eccentricities for the 
largest objects (e.g., Wetherill & Stewart 1993; Agnor & Ward 2002; Kominami & Ida 
2002 ). 

Figure 7 illustrates the evolution of the semimajor axes of the oligarchs for three cal¬ 
culations with So = 2 g cm“^ (lower panel), Sq = 4 g cm“^ (middle panel), and Sq = 8 g 
cm“^ (upper panel). As in Figure 1, the tracks change color when two oligarchs collide and 
merge. The labels indicate the final mass, in Earth masses, of the largest oligarchs at 100 
Myr. 

The timescale for the onset of chaotic growth depends on the initial surface density. 
More massive disks reach the transition first, (see, for example, Lissauer 1987). For Sq = 
8 g cm“^, the transition begins at ~ a few x 10® yr. The transition is delayed to ~ 1 Myr 
for Sq = 2 g cm“^. 

The character of the transition to chaotic growth also depends on the initial surface 
density. In relatively massive disks with Sq ~ 8 g cm“^, many oligarchs develop highly 
eccentric orbits and exhibit large variations in their semimajor axes. These large excursions 
result in many mergers and a rapid reduction in No. In less massive disks with Sq < 2-4 
g cm“^, only 1 or 2 oligarchs develop highly eccentric orbits. Most mergers are caused by 
two-body interactions, instead of large-scale dynamical interactions throughout the grid. 

Figures 8-10 illustrate these general conclusions. In Figure 8, the orbit crossing param¬ 
eter rapidly approaches zero for calculations with Sq = 8 g cm“^. At 0.1-1 Myr, po has a 
long plateau; close approaches between oligarchs cause po to fall below zero; mergers cause 
Po to jump above zero. During a series of 4 mergers at 10 Myr, po remains below 0 for a 
long period. After the hnal merger, Po jumps to 40, where it remains for many Myr. For 
smaller Sq, Po remains well above zero until one or two close pairwise interactions pushes it 
below zero. Once these interactions produce a merger, the systems stabilize and Po moves 
well above zero. 


Figure 9 shows the time evolution of the total number of n-bodies within the grid. 
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During the transition from runaway growth to oligarchic growth, No peaks and remains 
constant for 0.1-1 Myr. Once orbital interactions begin, No declines. The rate of decline is 
fastest in the most massive systems. 

The Hill parameter also shows the rapid transition from oligarchy to chaos. When No 
peaks, ph reaches a maximum. It declines sharply when mergers rednce the nnmber of 
oligarchs. After a merger or series of mergers, pn rises slowly as the remaining oligarchs 
accrete leftover planetesimals. 

To examine the sensitivity of these resnlts to rupro, we consider a broad range of pro¬ 
motion mass. For Sq = 1-16 g cm“^ at 0.8-1.2 AU, calculations with nipro ^ 2 x 10^^ 
(So/4 g cm“^) g cannot follow the evolution of the most massive coagnlation particles accn- 
rately. Thus, the production rate of oligarchs is inconsistent. Once a few oligarchs form, the 
models provide a poor treatment of interactions among oligarchs and therefore yield poor 
estimates of the merger rate and the timescale for the transition from oligarchy to chaos. 

Calculations with smaller nipro provide better solntions to the evolution of the largest 
objects. To measure the quality of the results for the transition from oligarchy to chaos, we 
dehne ts as the time when oligarchs with m > 10^®(So/8 g cm“^) g contain half of the initial 
mass and tm as the time when these oligarchs start to merge and the nnmber of oligarchs 
starts to decline. The scaling of the oligarch mass with Sq provides a clean way to compare 
models with different initial conditions. 

Our calculations show clear trends for tm and ts as a fnnction of Sq (Figure 11). The 
timescale for oligarchs to contain half of the initial mass ronghly varies as oc Tiq ' with 
very little dispersion^. The transition from oligarchic growth to chaotic grwoth oX t = tm 
also depends on Sq. Onr resnlts yield an approximate relation, tm oc However, this 

trend appears to have inflection points for S < 2 g cm“^ and S > 8 g cm“^. Calculations in 
progress will allow us to test this relation and its origin in more detail. 

There are no large trends for tm and fs with mpro- Although calculations with smaller 
nipro yield smaller dispersions in tm and t-£ at each Sq, the median valnes for tm and t^ are 
fairly independent of nipro- A Spearman rank test (Press et al. 1992) snggests a weak inverse 
correlation between ts and nipro, which we plan to test with additional calculations. 

Independent of nipro, onr results demonstrate a clear trend of the ratio = tm/ts with 


^The origin of this relation lies in the growth and stirring rates. In the absence of stirring, the growth 
time is t oc Sq Because more massive oligarchs form in more massive disks, the vertical scale height H of 
leftover planetesimals is larger in more massive disks. Larger scale heights reduce the growth rate by roughly 
oc which results in oc 
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initial surface density (Table 1). In massive disks with Sq ^ 12 g cm“^, mergers among 
oligarchs begin before oligarchs contain half of the mass ^1)- In low mass disks with 

^0 ^ 2 g cm“^, oligarchs start to merge after they contain half of the mass ^ !)• 

To test whether variations in the frequency and strength of dynamical interactions 
among oligarchs cause the trends in tm and r^, we dehne a ‘nearest neighbor parameter’ Un 
which measures the average number of oligarchs within 10 Rh of another oligarch. From 
§3.2, configurations with > 1 lead to strong dynamical interactions among the oligarchs; 
oligarchs interact mildly when < 1. Table 1 lists the average of the maximum value of 
Un and its dispersion for our calculations. Disks with > 0.5 have nn,max ^ 1; disks with 
^ 0.5 have Un^max ^ 1- This result conhrms the visual impressions from Figures 7-10: 
dynamical interactions among oligarchs are stronger in massive disks and milder in low mass 
disks. 

To conclude this section, Figure 12 shows the evolution of the mass for each oligarch 
in one calculation. Here, points of one color correspond to the track of one oligarch. From 
~ 10^ yr to ~ 10^ yr, large objects form and grow rapidly. Once stirring reduces gravitational 
focusing factors, growth slows. During the late stages of runaway growth and the early 
stages of oligarchic growth, several neighboring oligarchs merge. As their gravitational reach 
extends, these larger oligarchs grow faster and faster. Smaller oligarchs cannot compete for 
leftover planetesimals and grow slowly. 

During chaotic growth, the largest oligarchs merge with smaller oligarchs. As the small 
oligarchs are depleted, the merger rate slows. With highly eccentric orbits, a few small 
oligarchs last for 10-30 Myr before colliding with a large oligarch. After 100 Myr, all but 
one small oligarch have collided and merged with the two large planets that remain at the 
end of the evolution. 


3.4. Planet formation at 0.4—2 AU 

These calculations begin with 5 km planetesimals in a torus extending from 0.4 AU to 2 
AU. We divide this region into 40 annuli and seed each annulus with planetesimals in nearly 
circular and coplanar orbits (cq = 10“® and io = To provide a contrast with previous 

simulations, the planetesimals have surface density S = So(a/l AU)“^, with Sq = 2-16 g 
cm“^ at 1 AU. Calculations with the more standard S cx yield similar results. As 

in §3.3, we do not consider fragmentation. This torus is larger than the 0.5-1.5 AU region 
examined by Weidenschilling et al. (1997, see also Kominami & Ida 2002), but similar to the 
torus described by Chambers (2001). Chambers (2001) starts his calculations with lunar 



- 15 - 


mass objects, instead of planetesimals. Bromley & Kenyon (2006) compare results derived 
from our hybrid code for calculations starting with 1-10 km planetesimals or 100-200 lunar 
mass objects. When we start with smaller objects, our calculations produce fewer oligarchs 
on shorter timescales than the Chambers (2001) calculations. 

Growth in the larger torus is similar to that in the smaller torus. For Sq = 8 g cm“^, 
mergers produce a few large objects with radii of ~ 10 km at 0.4 AU in roughly a thousand 
years. Once runaway growth begins, it takes only ~ 10^ yr for the hrst promotion into 
the n-body code. As runaway growth propagates outward in heliocentric distance, objects 
throughout the grid reach the isolation mass of ~ 10^® g. It takes less than 10^ yr to produce 
10 isolated objects and another 2 x 10^ yr to produce the next 10 isolated objects. These 
objects continue to grow and reach typical masses of 0.05-0.1 in ~ 1 Myr. 

During oligarchic growth, occasional two body interactions produce mergers of oligarchs. 
These mergers always occur when the coagulation code produces 2-3 oligarchs in the same 
annulus. Often, dynamical interactions move one or two oligarchs into neighboring annuli, 
where they begin to accrete planetesimals more rapidly. Occasionally, dynamical interactions 
between several oligarchs in neighboring annuli produce one or two mergers, which stabilizes 
the system locally and leads to the isolation of the remaining 2 or 3 oligarchs. 

As oligarchs start to form at the outer edge of the grid, oligarchs at the inner edge of 
the grid begin to interact dynamically (Figure 12). For Sq = 8 g cm“^, these interactions 
start at ~ 1 Myr. Chaotic interactions then spread throughout the grid. For the next ~ 
10 Myr, mergers produce fewer but larger oligarchs, which rapidly sweep up the remaining 
planetesimals. After ~ 200 Myr, only 5 oligarchs remain. The largest have masses compara¬ 
ble to those of the Earth and Venus. Smaller oligarchs have masses comparable to the mass 
of Mars. 

As in §3.3, the evolution in less massive disks proceeds more slowly and less chaotically. 
The magnitude of dynamical orbital interactions is non-linear: once chaos starts in a massive 
disk, it rapidly propagates throughout the grid. Pairwise interactions dominate the dynamics 
of lower mass disks. These interactions usually do not affect other oligarchs. 

Figures 13-14 show the time evolution of po and pn- In all calculations, the typical 
orbital separation of the largest objects rapidly approaches zero when the hrst oligarchs 
form (Figure 13). This transition occurs sooner in more massive disks. During oligarchic 
growth, the orbits of oligarchs slowly move closer and closer. The number of oligarchs and 
the average Hill radius of the oligarchs (Figure 14) rise dramatically; these peak when 
~ 0.4-0.5. Once > 0.5, chaotic interactions cause oligarchs to merge. When the 

orbits overlap, large-scale chaos leads to rapid mergers and planets with masses comparable 
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to the mass of the Earth. In less massive disks with Sq < 2 g cm“^, the oligarchs remain 
well-separated relative to their Hill radii (Figure 14). Pairwise interactions between oligarchs 
then produce most of the mergers. This evolution is slow and results in more planets with 
masses comparable to the mass of Mercury. 

Finally, Figure 15 shows the growth of the oligarchs for a calculation with Sq = 8 g 
cm“^. Here, tracks of a single color show the time evolution in the mass of a single oligarch. 
Tracks end when an oligarch merges with another oligarch in the grid. 

The tracks in Figure 15 clearly illustrate the three stages of growth in the terrestrial 
zone. Most tracks are initially steep and then slowly turn over. The dramatic changes in 
the slope mark the transition from runaway growth to oligarchic growth. At 10® yr, this 
transition propagates as a wave from the inner part of the grid at 0.4 AU to the outer part 
of the grid at 2 AU. 

The onset of distinct steps in mass indicates the transition from oligarchic growth to 
chaotic growth. Starting at ~ 1~3 Myr, this transition tales ~ 3-10 Myr to propagate from 
0.4 AU to 2 AU. Once the full grid is chaotic, mergers rapidly separate the oligarchs into 
small oligarchs, with m ~ 0.01-0.1 ttiq and e > 0.1, and large oligarchs, with m > 0.2-1 me 
and e < 0.05. After ~ 100 Myr, this division is stark: there are 3 oligarchs with m > 0.5 m® 
and 3 others with m < 0.07 m®. 

Calculations with a large range in nipro provide a measure of the quality of these con¬ 
clusions. Because the 0.4-2 AU grid is much larger than the 0.84-1.16 AU grid, these 
calculations produce more oligarchs and are less sensitive to rupro than calculations in a 
smaller grid. Nevertheless, calculations with nipro ^ 3 x 10^® (So/4 g cm“^) g yield a much 
larger range in the evolution timescales than calculations with smaller nipro- As in §3.3, 
models with nipro ~ 1 — 10 x 10^^ g yield consistent results for these timescales. 


3.5. Other Physics 

In this set of calculations with the hybrid code, we concentrate on mergers and dynamical 
interactions between solid objects. We ignore fragmentation and interactions with the gas. 
In previous calculations, fragmentation of 1-10 km planetesimals produces cm- to m-sized 
fragments which are more closely coupled to the gas than larger planetesimals. Gas drag 
circularizes the orbits of these objects, which can then be accreted more rapidly than the 
leftover planetesimals. From previous numerical calculations, the onset of runaway and 
oligarchic growth are ~ 25% sooner than illustrated in §3.3 and §3.4 (e.g., Wetherill & 
Stewart 1993; Kenyon & Luu 1999; Kenyon & Bromley 2004b, 2005). 



Fragmentation is important for setting the visibility of the disk (Kenyon & Bromley 
2004b). The rate of debris production from fragmentation typically peaks during the tran¬ 
sition from runaway to oligarchic growth. If fragmentation produces an efficient collisional 
cascade, radiation pressure and Poynting-Robert son drag may eject small particles before 
the oligarchs can accrete the fragments. The disk then produces a substantial infrared excess 
(Kenyon & Bromley 2004b). If some process halts or slows the cascade, oligarchs can ac¬ 
crete the fragments efficiently (e.g., Goldreich, Lithwick & Sari 2004b). The disk may then 
produce a modest infrared excess. We have started calculations to evaluate the importance 
of fragmentation in setting the timescales for runaway and oligarchic growth and to estimate 
the amount of dusty debris as a function of time. 

Interactions with the gas can produce signihcant evolution in the eccentricity and orbital 
semimajor axis of oligarchs (e.g., Artymowicz 1993; Agnor & Ward 2002; Kominami & 
Ida 2002; Tanaka, Takeuchi & Ward 2002; Tanaka & Ward 2004). In the calculations 
here, simple gas drag as in Adachi et ah (1976) has a negligible impact on the evolution 
of planetesimals and oligarchs. However, interactions with density waves can circularize 
the orbits of oligarchs on short timescales (e.g., Artymowicz 1993; Agnor & Ward 2002; 
Kominami & Ida 2002; Tanaka & Ward 2004) and can cause signihcant inward migration 
of the orbit (e.g. Artymowicz 1993; Tanaka, Takeuchi & Ward 2002). For a gaseous disk 
with the surface density of a minimum mass solar nebula, the circularization timescale for 
an oligarch is Tc ~ 0.1-0.2 Myr at 1 AU (Agnor & Ward 2002; Tanaka & Ward 2004). 
Signihcant orbital migration of an oligarch can occur over ~ 1-3 Myr (Tanaka, Takeuchi 
& Ward 2002; Tanaka & Ward 2004). Thus, interactions with the gas occur on timescales 
comparable with timescales for oligarchs to grow and interact dynamically. 

The estimates for orbital migration and eccentricity damping assume a large gas density 
in the disk. Observations of young stars suggest the dust - and presumably the gas - 
disappears in ~ 1-10 Myr (e.g., Calvet et ah 2005; D’Alessio et ah 2005, and references 
therein). If the gas density declines exponentially with time as in equation (7) with to ~ 
1 Myr, radial migration may have little impact on the evolution of oligarchs. Eccentricity 
damping, however, may play an important role in the transition from oligarchic to chaotic 
growth. We have begun calculations to see how eccentricity damping and radial migration 
change the outcomes of our calculations. 

Together, fragmentation and interactions with the gas probably set the eccentricity of 
the hnal ensemble of planets. If fragmentation is efficient at converting intermediate mass 
objects into smaller fragments, the fragments can efficiently circularize the orbits of growing 
oligarchs (equation 2; Figures 5-6). Interactions with the gas also circularizes the orbits of 
the most massive objects. We have started a set of calculations to test the relative efficiency 



of damping and fragmentation in setting the eccentricity of surviving oligarchs. 


4. SUMMARY 

Our calculations with the hybrid n-body-coagulation code are the hrst to evolve an 
ensemble of planetesimals into a system of a few planets. The calculations reproduce both 
the standard results of coagulation calculations for the early stages of planet formation and 
the results of n-body calculations for the late stages of planet formation (see also Bromley 
& Kenyon 2006). In particular, we follow the evolution through the critical transition 
from oligarchic growth to chaotic growth, where objects evolve from isolated oligarchs into 
full-fledged planets. These calculations provide some new insights into the transition from 
oligarchic growth to chaotic growth (see also Kokubo & Ida 2002; Kominami & Ida 2002). 

In a disk with initial surface density Sq > 3-16 g cm“^ at 1 AU, collisions and mergers 
of 1-10 km planetesimals naturally lead to the formation of terrestrial planets with masses 
ranging from 0.05-2 (see also Wetherill & Stewart 1993; Weidenschilling et al. 1997; 
Chambers 2001; Kokubo & Ida 2002, and references therein). The growth follows a standard 
pattern. Slow orderly growth rapidly gives way to runaway growth, which concentrates 
much of the initial mass into many protoplanets with masses ~ 10^® g. Stirring of leftover 
planetesimals by the largest objects reduces gravitational focusing factors and slows growth. 
During oligarchic growth, large oligarchs become isolated and slowly accrete the leftovers. 

Several factors produce a transition from oligarchic growth to chaotic growth (see also 
Kominami & Ida 2002). As the oligarchs grow, they contain an ever increasing fraction of 
the total mass. A few oligarchs merge but most remain isolated from other oligarchs. As 
their Hill radii grow, their orbits begin to overlap. When the surface density in oligarchs 
reaches a critical value, oligarchs interact chaotically. As gravitational interactions scatter 
oligarchs throughout the disk, the merger rate increases dramatically. Eventually only a few 
oligarchs remain in roughly circular orbits. 

Our results isolate the two conditions necessary for the transition from oligarchy to 
chaos. When oligarchs contain roughly half of the mass of solid material in the disk, 
dynamical interactions between oligarchs are more important than dynamical friction from 
planetesimals (e.g., Goldreich, Lithwick & Sari 2004b). When the surface density in oligarchs 
exceeds a critical value, ~ 2-3 g cm“^, oligarchs begin to interact chaotically (§3). More 
massive disks can reach this limit when S/ < (Table 1). In less massive disks, milder 
dynamical interactions begin when (Figure 11). If the surface density in oligarchs 

remains below the critical value, interactions between oligarchs are less chaotic even when 
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S; > Ss- Interactions among 2 or 3 oligarchs produce a small merger rate which eventually 
yields a system with lower mass planets compared to more massive disks. 

Although dynamical interactions among the ensemble of oligarchs produce terrestrial 
mass planets in all disks, more massive disks yield more massive planets. Our results suggest 
a maximum mass, nimax ~ 1“2 m® for Sq ~ 8-16 g cm“^ and nimax ~ 0.1-0.2 m® for Sq ~ 
1-2 g cm“^. Because young stars appear to have a wide range of initial disk masses, we expect 
a wide range in the masses of terrestrial planets in exosolar systems. In future studies, we 
plan to address this issue in more detail. 

In terrestrial planet formation, the transitions between different stages of growth pro¬ 
duce distinct waves through the disk. During the transition from orderly growth to runaway 
growth, the increase in the collision rate rapidly propagates from the inner edge of the disk 
to the outer edge. This transition is rapid and takes < 10® yr to move from 0.4 AU to 2 AU. 

Although less rapid, the transitions from runaway to oligarchic growth and from oli¬ 
garchic to chaotic growth also propagate from the inner disk to the outer disk. During the 
transition to chaotic growth, dynamical interactions tend to produce more chaotic orbits at 
the outer edge of the disk. This behavior depends on the surface density gradient. In disks 
with steep density gradients, S ~ Soa“"' with n >3/2, chaotic interactions propagate slowly 
outward. In disks with shallower density gradients, n < 1, dynamical interactions tend to 
concentrate more mass in the outer part of the disk. This difference in behavior is set by the 
growth rate, ~ P/S ~ where P is the orbital period: the wave of growth propagates 

more rapidly through disks with shallower density gradients (e.g. Lissauer 1987). 

These results have several interesting consequences for the evolution of planets in the 
terrestrial zone (see also Kominami & Ida 2002). The transition from oligarchic to chaotic 
growth occurs on timescales, ~ a few Myr, well before radiometric evidence suggests the 
formation of the Earth was fairly complete (e.g., Yin et al. 2002). Planets are also fully 
formed well before the estimated time of the Late Heavy Bombardment, ~ 100-300 Myr 
after the formation of the Sun (e.g., Tera, Papanastassiou, & Wasserburg 1974; Hartmann 
1980; Ryder 2002; Koeberl 2003) 

Throughout the chaotic growth phase, our calculations produce many lunar- to Mars¬ 
sized objects on highly eccentric orbits. These objects are good candidates for the ‘giant 
impactor’ that collided with the Earth to produce the Moon (Hartmann & Davis 1975; 
Cameron & Ward 1976; Benz, Slattery, & Cameron 1986; Canup 2004a,b). As we complete 
calculations with fragmentation and migration, predicted mass and eccentricity distributions 
for these objects will yield better estimates for the probability of these events. 

Together with the dynamical influence of Jupiter (e.g., Kominami & Ida 2004), the 
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highly eccentric orbits of lower mass oligarchs have an important role in clearing the inner 
solar system (e.g., Goldreich, Lithwick & Sari 2004a; Kominami & Ida 2004) and delivering 
water to the Earth (Lunine et ah 2003). Traditionally, data from our solar system provide 
the only tests of clearing mechanisms (e.g., Grogan et ah 2001; Nesvorny et al. 2002a,b, and 
references therein). In the next decade, comparisons between predicted infrared excesses from 
the debris disks leftover from terrestrial planet formation and observations from Spitzer and 
TPF-Darwin will yield new constraints on clearing timescales and the evolution of volatile 
species in the terrestrial zone (e.g., Beichman et al. 2005). These comparisons will enable 
better numerical calculations and an improved understanding of the hnal stages of terrestrial 
planet formation. 
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Table 1. Transition from Oligarchy to Chaos at 0.84-1.16 AU 


So (g cm ^ 

0 log (yr) 

tt. 

^n,max 

1 

7.0 ± 0.05 

0.58 ± 0.01 

0.85 ± 0.05 

2 

6.6 ± 0.05 

0.53 ± 0.01 

0.91 ± 0.05 

4 

6.3 ± 0.07 

0.48 ± 0.02 

0.99 ± 0.10 

8 

6.0 ± 0.08 

0.41 ± 0.02 

1.11 ± 0.07 

12 

5.7 ± 0.10 

0.32 ± 0.03 

1.24 ± 0.12 

16 

5.3 ± 0.10 

0.24 ± 0.03 

1.33 ± 0.11 
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Fig. 1.— Evolution of eccentricity for oligarchs in a swarm of lower mass planetesimals. 
Each colored track shows the eccentricity of a single oligarch with time. A change in the 
color of a track (for example, in the middle panel at 1 Myr when the dark gray track turns 
blue) indicates a merger of one pair of oligarchs. The legend of each panel lists the surface 
density in g cm“^ at 1 AU of oligarchs with individual masses rus = 10^® g. In each panel, 
the surface density of planetesimals with masses m<j = 10^^ g equals the surface density 
of oligarchs. When < 1 g cm“^, the orbits of oligarchs do not overlap on timescales of 
1 Myr. For S; > 2 g cm“^, gravitational stirring leads to overlapping orbits and mergers 
of oligarchs. More massive systems develop overlapping orbits on shorter timescales, which 
leads to more mergers on timescales of 1 Myr or less. 
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Fig. 2.— As in Figure 1 for the semimajor axes of the oligarchs. The legend in the lower 
left corner of each panel indicates the number of oligarchs at the beginning and end of each 
calculation. 
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Fig. 4.— Evolution of the orbit crossing parameter po for models with = 1-8 g cm“^. 
For Po ^ ‘2^, orbits of oligarchs do not overlap and mergers do not occur. As po approaches 
0, overlapping orbits of larger objects lead to mergers. 
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Fig. 5.— Evolution of semimajor axes for oligarchs in a swarm of lower mass planetesimals. 
As in Figure 1, each colored track show the position of one oligarch; changes in color along 
a track indicate a merger. The legend of each panel lists the number of oligarchs and the 
surface densities of planetesimals with = 10^^ g (S^) and oligarchs with mi = 10^® g (S;) 
at 1 AU. As more mass is placed in smaller objects, dynamical friction Erst reduces orbit 
overlap between oligarchs. However, when the planetesimals contain more mass, viscous 
stirring between planetesimals aids orbit overlap and leads to mergers of the oligarchs. 
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Fig. 6.— As in Figure 4 for = 10^® g. Viscous stirring between planetesimals is minimal 
and does not lead to overlapping orbits of oligarchs in 1 Myr. 
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Fig. 7.— Evolution of semimajor axes for oligarchs with masses larger than rripro = 
10^® (Eo/4 g cm“^) g in a full planet formation calculation at 0.86-1.14 AU. The legend 
of each panel lists the initial surface density in 1-10 km planetesimals at 1 AU (upper left 
corner, in units of g cm“^), and the mass of the largest planets (in m®; at the right end of 
each track). The dashed lines indicate the extent of the planetesimal grid in the coagulation 
code. Runaway growth produces protoplanets in 10^ yr to 10^ yr. Continued growth of 
protoplanets during the oligarchic phase leads to orbit overlap on timescales of roughly 1 
Myr. More massive systems reach this limit faster than less massive systems. Overlapping 
orbits leads to merger and more rapid growth of protoplanets. 
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Fig. 8.— Evolution of the orbit crossing parameter po for the models of Figure 7. The 
legend lists the surface density in units of g cm“^ for each model. The dashed lines indicate 
when the surface density of oligarchs is half the total surface density. The leftmost dashed 
line corresponds to the largest initial surface density; the rightmost line corresponds to the 
smallest initial surface density. As oligarchs form, po declines and reaches a plateau where 
Po ~ 1-5. Continued growth of the oligarchs leads to a minimum in po, roughly when 
Tii ~ Tis- minima Subsequent mergers reduce po until the system is in rough equilibrium, 
when Po increases dramatically. 
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Fig. 9.— Evolution of the number of oligarchs No for the models of Figure 7. 
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Fig. 10.— Evolution of the Hill parameter ph for the models of Figure 7. The dashed lines 
indicate when the surface density of oligarchs is half the total surface density. The leftmost 
dashed line corresponds to the largest initial surface density; the rightmost line corresponds 
to the smallest initial surface density. When the number of oligarchs is maximum, the Hill 
parameter peaks. Mergers reduce pn] growth increases pn- Chaotic growth occurs when 
PH > 0.07-0.1. 
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Fig. 11.— Variation of and tm as a function of Sq for different values of nipro- The legend 
indicates symbols for the logarithm of the promotion mass. For each promotion mass and 
surface density, the points indicate typical values and extreme values. Thus, the two panels 
provide a visual impression of the dispersion of the results among a large set of calculations. 
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logE(g cm'^) 

Fig. 12.— Evolution of masses for oligarchs with masses larger than 2 x 10^^ g in a full planet 
formation calculation at 0.86-1.14 AU. The legend lists the initial surface density in 1-10 
km planetesimals at 1 AU (upper left corner, in units of g cm“^). Each colored track shows 
the mass evolution for one oligarch. Discontinuities or terminations in the tracks indicate 
mergers of large objects. 
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— As in Figure 7 for a calculation at 0.4-2 AU with rupro = 2 x 10^®(S/4 g cm g. 


Fig. 13. 
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Fig. 14.— As in Figure 8 for the models of Figure 13. 
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Fig. 15.— As in Figure 10 for the models of Figure 13. 
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Fig. 16.— As in Figure 12 for a model from Figure 13. 



